Turbulent behaviour in magnetic hydrodynamics is not universal 
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■ Abstract 
O . 

^ , A short distance expansion method (SDE) that is weU known in the quantum 

field theory for analysis of turbulent behaviour of stochastic magnetic 
^ hydrodynamics of incompressible conductive fluid is applied. As a result is shown 

■ that in an inertial range the turbulent spectra of magnetic hydrodynamics depend 
Q ! on a scale of arising of curls. 

> 

H ■ 1 Introduction 

, 

The methods of quantum fleld theory are successfully applied for description 
of the critical phenomena and developed turbulence recently. This approach has 
an important advantage before the classical one. For example it enables one to 
renormalize all the correlation functions of a model as well as to deflne their various 
asymptotics. 

A stochastic problem of magnetic hydrodynamics (MHD) of incompressible 
conductive fluid with external random force and rotor of random current 

VtVi = uAvi - d,p + {63)6, + Vt = dt + vd, (1) 

VA = y'^6, + {6^)v^ + f! 

is equivalent to a theory of four flelds with generating functional of renormalized 
correlation functions (Green functions) 

G{A^) = J V^detMexp{SR{^) + ^A^), $ = {^,^',6,6'} 
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in which the renormahzed functional of action is 

+^'[-dt^ + Ziz/A^ - {^d)^ + z^{ed)e] + e'[-dte + Z2vuM - {^d)e + {eo)^] 

(necessary integration on x and t as well as summation on repeated badges are 
implied). Here the fields (p{'x,t) and 6'(x, t) are both transversal (the vector field of 
velocity ip due to incompressibility of fluid {di(pi = 0) and the pseudovector field 9 
as connected with transversal field of a magnetic induction B : = B/^/47Fg {g is 
a density of fluid, z/ is a viscosity, p is a pressure). We shall use the dimensionless 
constant A (it means inverted number of Prandtl) instead of u' = c^/Ana ( here a 
is a conductivity, c is the velocity of light) by this way u' = Xu. 

Under construction of the model existence of an inertial range is implied therein 
a real external energy pumping (correlators of the random forces) can be simulated 
by exponential model of (5-function. 

= QlO^oPisdwi D^is = 920^0 -Pisd-ee , -Djf = 930^0 ^ism^md-ve, (3) 

d^^ = k'-'-'^, dee = k'-'-'"' , d^e = k'-"-^'^^^^ . 

In momentum-frequency representation (from frequency the correlators Di^ do 
not depend) Pis is the transversal projector, d means the dimension of space (com- 
pletely antisymmetric pseudotensor Sism is determined only at d = 3). 

The factors go in correlators are played a role of charges; the positive constant a 
is an arbitrary parameter. The parameter e serves for construction of 
decomposition of correlation functions, as the physical value of e the = 2 is 
considered that simulates energy pumping from the large-scale movements of fluid. 

The model (|I] was investigated in in the first order on e with the help of 
recursive renormalization group method. It is shown that in the system two different 
regimes of developed turbulence can be realized: "kinetic" and "magnetic" that are 
connected with existence of two infrared (IR) stable fixed points of RG transforma- 
tion. In the critical asymptotics of spectral density of energy < y9(k)(y9(— k) > 
and < ^(k)6'(— k) > were determined. In more general formulation of the prob- 
lem with inclusion of the cross correlator of random force D*^^ was considered by 
the quantum field RG method. The renormalization constants Zi were calculated 
there in the first order of (^q,: 

^ _ ^ gid{d-l) g2id^ + d-A) ^ _ ^ Qi 92 
^ 4:Be ABaX^e ' ^ BXe BaX^e' 

gi{d + 2){d-l) _ g,{d + 2){d-3) _ 
^ 2BX{X + l)e 25aA2(A + l)e' 

here B = d{d + 2)(47r)'^^^r((i/2) and r(a;) is the gamma-function. 
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Renormalizability of the model was proven, and situation in the charging space 
of IR-stable fixed points corresponding to kinetic and magnetic regimes were found: 



' - 3u e{l + X^) \/43/3 - 1 
^?u = ^ = ^^, 92. = 0, K = , (4) 

( the region of stabihty of the point is a < 1.16. ) 

9u = A, = 0, = = «^ (5) 

( this one is stable at a > 0.25). In paper features of a scaling behaviour in 
the model were investigated, and critical dimensions of fields and parameters of the 
theory are found out in the both critical regimes. 

We shall be interested in dependence of characteristics of developed turbulence 
(the correlation functions) in the inertial range from conditions of arising of large- 
scale curls. In the theory of developed turbulence it is supposed, that the energy 
pumping into the inertial range is executed by the vortices of a large size A. We 
shall take into account this scale having supplied the model (|^ ) by an infrared 
mass parameter m = 1/A to consider the relative corrections to developed turbulent 
spectra with the help of SDE method. 

According to this method following the operator decomposition is fair: 

0(xi, t)0(x2, t) ^ 5] Q(r)F,(x, t); (6) 

i 

where x = (xi + X2)/2, r = xi — X2, Fi are various local averages (composite 
operators) . The averaging of (^ yields the asymptotics for pair correlation functions 
at mr 

< 0(xi,t)0(x2,t) >~ C^) 

i 

here Ap. are the critical dimensions of the composite operators; means some 
constants. On decomposition (|^) we conclude, that from the point of view of an 
opportunity of transition to a massless theory ( m — ) in the inertial range the 
operators with a negative critical dimension are dangerous. 



2 Critical dimensions of the composite operators in the model 
of magnetic hydrodynamics 

The renormalized operators are defined by the formula Fi = ZikF^{<^pt). On the 
known matrix Zik a matrix of anomalous dimensions 7j/c = {Z~^)ijVMZjk is calcu- 
lated and then the matrix of critical dimensions Aik = {dp)ik + Ai^{dp)ik + 'yik where 
Ai_j designates critical dimension of frequency, and dp, d% are the momentum and 
frequency canonical dimensions of F; they are being determined from requirement 
of frequency and momentum dimensionless of terms of the action functional. 



3 



The particular critical dimensions are eigenvalues of the matrix Aj^. They corre- 
spond to linear combinations of composite operators Lj(F^) = UikF^ which diago- 
nalize the matrix Ajfc. 

We shall consider the dimensions of elementary composite operators of the MHD 
model: 0i0A:, ^is well as vector operators {dcjxp)^ with various transpositions of 
badges. 

The tensor is a sum of two independent tensors 

Convolution of the first operator on badges yields family of scalar operators 

^1 = ]^ViVu F2 = ^OiOf, F3 = OiVi. 

Trace of the second expression in (H) is equal to zero. 

Essential property of the theory (|I| - |^) is Galileian invariancy; it was used in 
m for investigation of the composite operators in a problem of usual stochastic 
hydrodynamics. A non-stationary Galileian transformations of the fields 

</?a(x, t) = (^(x + u{t),t) - a(t); (^;,(x, t) = v?'(x + u{t),t); 

Oai^, t) = ^(x + u(t), t); ^^;(x, t) = e'{^ + u(t), t); 

(a parameter of transformation - v(t) is the vector function dependent only from 
a time well decreasing at \t\ 00 and u(t) = dtw(t') ) reahzed in the Ward 
identity leads [^j to 

Jdx{ao.[^^ - g^^d^Vk] + dAao.g^^]} < 00. (9) 

Here aoa are the nonrenormalized functions of sources of family of operators Fa, and 
symbol "< cxd " means that the considered functional is finite. 

We shall assign Gi = ViVk/2 and G2 = 9^9^/2. The operator G3 = Vi9k is pseu- 
dotensor so it doesn't mix with the first two operators due to renormalization. Gi 
is noninvariant of the Galileian transformations so the formula permits to prove 
that Gi is finite {Zu = 1), and that it doesn't mix to the operator G2 in proceeding 
of renormalization, hence Z21 = (one can say, that the operator Gi aren't being 
renormalized and doesn't mix to the Galileian invariant operator G2 since Gi is not 
Galileian invariant). Besides, from (j^) the absence of divergences in diagrams with 
the composite operator G3 (it is noninvariant of Galileian transformations also) is 
followed. Then from the definition of renormalized composite operators it is easy to 
get the expression: Z33 = Z^^g^ = Zg^. 

Thus, the matrix of renormalization constants Z^s of the operators Gi is 

= 



1 


Zl2 








Z22 
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Unknown elements of the matrix one can calculate with the help of standard dia- 
grams with the following propagators 

< <^,(k,t)^;(-k,0) >o= p^.e-^^^'eit); 
< Viik,t)^,{-k,0) >o= ^P,,VA;'-"-'^M2-e-'^'='l*l; 

< ^,(k,t)^;.(-k,0) >o= P,,9{t)e-''^'"'- 

1 „t 2l.2-d-2ae 

< 9,ik,t)e,i-k,0) >o= -P^,^ M^-e-^^'^ 1*1; (10) 

<9,{k,t)^,{-k,0) >o= g.s /''''^' ' ^'^"^' M(^+")^fc,[e-"-^"*^(t) +e-^'*g(-t)]; 

[i + A) 

Account of the renormalization constants for the scalar operator we execute having 
curtailed the badges 6ik. 

The calculations in one-loop approximation under diagrams specified on a fig. 1 
of the appropriate constants Z^s in the theory (H) yelds a renormalized action for 
the generating functional of correlation functions with the composite operators F^. 
Sr = Sr + aoaZasFs^{Z^(j)); here the renormalization constants are 



2(A + 1) ^7~^s^' ^ - ^ + 2(A + 1) ^7 " '^e' 



Z'^^'^12 = - - -), Z'^'^22Zl = 1 + - ^).(ll) 

2(A+1)^£ ae^' ^ 2(A + l)^e ae' ^ ' 

As far as in the fixed points of RG transformation some values of the charges 
approach to zero in avoidance of trivialization of asymptotics in all orders of the 
^-decomposition it is useful to redefine the fields as follows: 



Q sJg2v'^M'''Q; V Jgiiy^'v; (12) 



The transformations don't change positions of the fixed points and values of the 
renormalization constants (the canonical dimensions of the operators Fa vary only). 

The critical dimensions of the operators Gi calculated on the constants Zaa in 
the first order of e and at any d are listed in table I. 

It should be noticed that the value of Ai is exact. The constant Z12 defines an 
admixture to Gi of the operator 6*2 • Considering (^) in a condition of the kinetic 



regime with the fields of (|T2D we have G2 — * 5'2G^2 and taking into account that in 
the kinetic fixed point (§) = i^'s easy to show that just the operators Gi, G2 
(instead of their mixture) have the dimension Aj in this point. 

The constant Z12 ~ 0(A) which is responsible for mixing of the operators disap- 
pears in the magnetic point. 
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We shall consider elements of a matrix Zik correspond to renormalization of the 
set of scalar and vector composite operators 0'0 and {d(j)(f))i. 

At this set there are the operators reducing to a total differential <9(00) with 
various transpositions of badges. Renormalization of them is equivalent to renor- 
malization of the operators 00 have considered above. The critical dimensions is 
being appropriated to these operators surpass dimensions located in table 1 per unit 
of. 

The operators of a 0'0-type don't mix to any other operators and aren't being 
renormalized because of 1-irreducible diagrams that is responsible for mixing of these 
operators with the other are equal to zero as far as they contain cycles of advancing 
lines. A structure of interactions in the (0) provides removal from each diagram 
a one derivative on each external line of a 0'-type that effectively lowers an index 
of divergence of the diagrams. Therefore, in the minimal subtraction (MS) scheme 
of renormalization that is appropriate to these operators the diagonal elements are 
Zaa = 1, and all nondiagonal ones are equal to zero. Thus, Afj,/^ = 3. 

The remaining vector operators Fi = Viv"^ and F2 = v6^ are true tensors, and 
-^3 = 6'j6'^, F4 = 9v'^ are pseudotensors; these pairs of the operators are being 
renormalized independently from each other. The operators -Fi_2 are noninvariant 
to Galileian transformations, thereof, it is easy to approve the finiteness of the 
operator as well as that it doesn't mix to vO"^ due to renormalization. Thus, 
Zn = I, Z21 = 0. 

Similarly, the Galileian invariant operator F3 can't mix to the operator F4 that 
means Z34 = 0. 

For definition of remaining elements of the matrix in a one-loop 
approximation it is necessary to consider the diagrams that is shown on a fig. 2. It 
gives 

_ ^ X{d+2){d-l) ^g[ g' (^+2)(rf^ g' 

^''--^ — — ^7"^^' (A+i) ^7'^^' 



(A 

\{d + 2){d-l) ,g[ g' 7 g[^g'2 

^43^0 — 7T — —r\ 1 j-i ^44^6* — J- 1 

[X + i)eae e ae 



The matrix of renormalization constants of this family of the operators has a block 

triangular form, so the critical dimensions are determined by the diagonal elements 
7 

The values of critical dimensions calculated on (pTSj) are listed in table 2. 

The nondiagonal elements Z12 and Z43 define an admixture of the operators Fi 
and F4 to F2 and F3. Taking into account that in the kinetic mode Fi = gi'^'^v^, 
F2 = gi'^g2v9'^ F3 = g2^^'^0^, and F4 = gig2^^'^dv'^ and that 5(2* = in the fixed 
point (IS) it is possible to assert the certain critical dimensions are belonged to the 
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operators, instead of their linear combinations. In the magnetic point (Q) Z\2 = 
Zm = 0. 



3 Discussion of the results 

It is important to note that any of the operators considered can't participate as 
amendments for phenomenological equations of MHD. The operators are being pos- 
sessed of the essential critical dimensions don't contain auxiliary fields 0', but the 
operators of a (j)'(f)-tjpe (a function of response) are inessential and don't satisfy the 
requirements of Galileian invariancy also. 

The operators of canonical dimension d = 3 define new nonanalytic corrections 
to the spectra of developed turbulence that was found in 0. According to SDE 
method such corrections for a pair correlation function < 0i02 > can be represented 
as follows: 

(here A,^ are the critical dimensions of the fields). 

As far as in the inertial range the following estimation is correct m/k <C 1, the 
formula (|T^ results to nonanalyticities in a case of A^. < 0. At the real value ofe = 2 
the critical dimensions of G2 (a > 1/2, Fi, F2 (a > 0.243), F3, F4 {a > 
0.82) in the kinetic point become negative. The dimension of F4 at a > 3/4 is 
negative in the magnetic point also. 

Thanking to Galileian invariancy of the theory we can refer to the results of 
where the terms of the sum of (0) for static correlation functions connected with the 
Galileian noninvariant composite operators were proven to yield not a contribution. 
Those in our case are Fi, F2, and F4. They can participate in the decomposition 
([T3|) only for dynamic correlation functions. The Galileian invariant operator F3 
gives the contribution in the sum in all the cases. 

Comparing of the dimensions of the composite operators of families 00 and (900) j 
at the real value e = 2 one can see that the set {d(f)(f))i appears more essential in the 
scaling range (so in the kinetic regime Aq2 > A^02 , Aga , A„2 > A^s ) . One can 
assume that the tendency of growth of essentiality of operators is being demonstrated 
by the elementary operators in the model of magnetic hydrodynamics will be saved 
for more complex operators. 

The results received by us testify that in the inertial range in MHD the correlation 
functions depend on the external integral turbulent scale m. Thus, the behaviour in 
the model is not universal. 

This important result can be checked experimentally besides the particular values 
of critical dimensions of the composite operators calculated here can be measured 
on an experiment too that will be served certainly to check of the offered theory. 
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Table I 



The critical dimensions of the composite operators of (f)i(f)j type 


Operator 


Kinetic point 


Magnetic point 




2 - 4/3e 


2 




2-2(a-3/10)£ 


2 + 3ae 




2- (a + l/3)£ 


2 + ae 


{GijijSij 


2 - 4/3£ 


2 


{G2)ijSij 


2 - 2ae 


2 + 12ae 


{GsjijSij 


2- (a + l/3)£ 


2 + ae 
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Table II 



The critical dimensions of the composite operators of {d(f)(f))i type 


Operator 


Kinetic point 


Magnetic point 




3-2e 


3 


F2 


3-2(a + 0.507)£ 


3 + 18ae 


F3 


3 - 3(a+ l)e 


3 + 60ae 


F4 


3- (a + 0.68)£ 


3-2ae 
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fig. 1 The set of one-loop 1-irreducible diagrams are responsible for the nontrival 
renormalization Constances of the composite operators of a 0^0^ type. 
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fig. 2 The set of one-loop 1-irreducible diagrams are responsible for the nontrival 
renormalization Constances of the composite operators of a type. 
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